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Abstract 

We consider a mass-less manifestly covariant linear Schrodinger equation. First, we show that it 
possesses a class of non-dispersive soliton solution with finite-size spatio-temporal support inside 
which the quantum amplitude satisfies the Klein-Gordon equation with finite emergent mass. We 
then proceed to interpret the soliton wave function as describing a particle with finite mass, energy 
and momentum. Inside the spatio-temporal support, the wave function shows spatio-temporal 
internal vibration with angular frequency and wave number that are determined by the energy- 
momentum of the particle as firstly conjectured by de Broglie. Imposing resonance of the inter- 
nal vibration inside the spatio-temporal support leads to Planck-Einstein quantization of energy- 
momentum. The first resonance mode is shown to recover the classical energy-momentum relation 
developed in special relativity. We further show that the linearity of the Schrodinger equation 
allows one to construct many solitons solution through superposition, each describing a particle 
with various masses, energies and momenta. 

PACS numbers: 03.65.Pm; 03.65.Ge 
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I. INTRODUCTION: RELATIVISTIC MADELUNG FLUID 



Surprisingly, despite of the remarkable pragmatical successes of quantum theory, there is 
still an unsettled debate concerning the physical status of its main ingredient, that is the 
wave function. Roughly speaking, there are two main attitudes toward this foundational 
issue flflSQl. 

The first cult considers that the wave function has no physical reality at all. They 
assume the wave function as the representation of our knowledge about the physical reality 
rather than to refer to the reality itself. This interpretation then regards quantum theory 
as a theory concerning our knowledge about the reality obtained through experiment rather 
than a theory about physical reality. It is developed by Bohr and Heisenberg and commonly 
known as Copenhagen interpretation. This line of thought eventually led to the probabilistic 
view of wave function through Born's rule j^]. 

The second attitude is to consider the wave function as a real physical field referring 
directly to the physical object being described, like say electromagnetic field. There are many 
interpretation of quantum theory which attribute such a physical status to the wave function. 
The mostly mentioned interpretations which support this view includes the axiomatic-mos ;- 
"used" standard interpretation of Dirac-von Neumann 6, 71, the Bohmian mechanics 4, 8|, 

un 

many worlds interpretation [9|, HQJ, the theory of spontaneous localization 



11], etc. 



The less mentioned one is de Broglie's theory of double solutions [12J, |_13(] . In his attempt 
to solve the dual nature of matter as particle and wave, he was searching for a nonlinear 
wave equation which assumes a non-dispersive soliton solution at the amplitude which is 
sufficiently high, while possesses a linear solution satisfying a linear superposition at the 
amplitude which is weak. He then proposed that the soliton part should be regarded as a 
particle which is guided by the linear part of the solution. This idea eventually led him to 
derive his famous guiding principle: 

E = hu, p = Kk, (1) 

which relates the energy-momentum {E, p} of the particle with the angular frequency-wave 
number {u, k} of the linear wave. 

Equations (0Q) can be argued as the most important principle of quantum mechanics [jjj . 
Partly inspired by those relations, Schrodinger developed his celebrated equation. Yet in 
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contrast to de Broglie's original idea, Schrodinger equation is linear with respect to the 
wave function. In this theory, a free particle is then usually represented by a plane wave 
satisfying the relations of Eqs. (pQ). Despite physically vague, the plane wave representation 
surprisingly works for all pragmatical purposes [3J. Yet one can ar gue that the successes of 



this representation relies heavily on the above de Broglie's relation 13) . It is apparently the 



pragmatical successes of the linear Schrodinger equation if combined with the Born's rule 



that eventually discourage people from further continuing de Broglie's program 



IS 



13, 



l^ | . In this paper, by considering a mass-less relativistic linear Schrodinger equation, 
we shall show that it has a new class of soliton solutions with properties exactly envisioned 
long time ago by de Broglie. 

Let us consider a closed system whose state is uniquely determined by a complex-valued 
wave function in spacetime, if)(q), where q = (q°, q 1 , q 2 , q 3 ) = (ct,x,y,z). Further, let us as- 
sume that the wave function satisfies the following manifestly covariant mass-less Schrodinger 
equation 

d ft 2 

ih—if)(q;X)^-Dif)(q;X). (2) 
Here, A is some affine parameter, □ = —r] ab d a db and r] ab = diag(— 1, 1, 1, 1) are 



D'Alembertian operator and flat Minkowskian metric, respectively. Eq. ([2 



posed to interpret the Klein-Gordon equation through a particle model 



2a 



has been pro- 



2l|. 



Next, putting the wave function into polar form, if) = I exp(iS/h), where I and S are 



real- valued functions, and separating into the real and imaginary parts, one obtains 



14: 



% = | +«.(,.•) = o. (3) 

Here, p(q; A) = \ip\ 2 = I 2 is the quantum probability density, v a (q; A) is a velocity field 
generated by the quantum phase S(q; A) as 

v a (q;\) = d a S(q;\), (4) 

and U(q; A) is generated by the quantum amplitude I(q; A) as 

UM - ^. (5, 



We have thus adopted the Madelung fluid picture for the Schrodinger equation 22]. Due to 



its formal similarity with the Euler equation in hydrodynamics, the term on the right hand 
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side of the left equation in Eqs. d3J), F a = —d a U, is called as quantum force field. Thus, 
correspondingly, U(q) is called as quantum potential. 

Let us remark that written in the form of Madelung fluid, it becomes clear that the original 
Schrodinger equation possesses a hidden self-referential property. Namely, the quantum 
potential U(q) is generated by the quantum probability density p(q) through Eq. dSJ). This 
in turn will dictate the way p(q) must evolve with time through Eqs. d2J) and so on and 
so forth. It is thus reasonable to expect some interesting self- organized physically relevant 
phenomena. In particular, in this paper we shall be interested to study the fixed points of 
such dynamics. 



II. SELF-TRAPPED QUANTUM PROBABILITY DENSITY 

Let us proceed to specify a class of wave functions whose quantum probability density is 



further related to its own quantum potential as 



23| : 



I f U(q) \ 



p(?) = Z(T) eXP l-^> (6) 
where T is a real- valued parameter below chosen to be non- negative and Z(T) is a normal- 
ization factor. We shall show that the above class of quantum probability densities possesses 
non-trivial and physically interesting properties. To do this, notice that combined with the 
definition of quantum potential given in Eq. (jSJ), Eq. (jSJ) comprises a differential equation 
for U(q) or p(q) subjected to the condition that p(q) must be normalized. In term of U(q), 



one has to solve the following nonlinear differential equation 



23| : 



~DU = ^d a Ud a U + ^U. (7) 

One observes that the above differential equation is invariant under Lorentz transformation. 
Hence, given a solution U(q), then any function U(q'), where q' a = A a b q b and A a b is Lorentz 
transformation, is also a solution of Eq. (jTJ). 

Let us develop a class of solutions in which the quantum probability density is being 



trapped by the quantum potential it itself generates |23|. To do this, let us assume that 
there is an inertial frame so that the quantum probability density is separable into its spatial 
and temporal parts as follows: 

P(q) =P*(x)Pt(t), x = {x,y,z}. (8) 



In this case, the quantum potential can then be decomposed into 



U(q) = U K (x) + U t (t) 1 (9) 

where 

Here, U = ^Jpi with i — x, i; and d 2 = <9 X ■ <9 X where <9 X = {c^, 9 y , d z }. 

The condition of Eq. (jSJ) is not Lorentz invariant so is the resulting class of solutions 
we are going to develop. Yet, its nontrivial property will be shown to be Lorentz invariant. 
Inserting Eq. ©, Eq. © can thus be re-collected as 9 x [/ x -(l/2T)9 x f/ x -9 x t/ x -(4T '/h 2 )U x = 
{l/c 2 )d 2 U t - (l/2T)(d t U t ) 2 + (4T/h 2 )U t = D, where D is constant. Below for simplicity we 
shall take the case when D = 0. One thus has to solve the following decoupled pair of 
nonlinear differential equations: 

1 4T 
d'iU x - ^d x U x ■ d x U x - — U x = 0, 



1 1 AT 

- 2 d 2 tU t --{d t U t ) 2 + -U t = V. (11) 



A. Spatial self-trapping 

Let us first discuss the spatial part by solving the upper differential equation in Eqs. (jTTT) . 
To do this, let us search for a class of solutions in which the spatial quantum probability 
density is further separable as 

Px(x) = p x {x)p y {y)p z {z), (12) 

so that the spatial part of the quantum potential is further decomposable into 

C/ X (x) = U x (x) + U y (y) + U z (z), Ui(i) = - y^, (13) 

where i — x,y, z. Putting this anzatz into the upper differential equation in Eqs. pip , one 
can choose a class of solutions in which each Ui(i) satisfies the following decoupled nonlinear 
differential equations: 

1 4T 

dfU i = —(d i U i ) 2 + —U i , i = x,y,z. (14) 

To avoid complicated notation, below we shall consider only the x— degree of freedom. The 
other two spatial degrees of freedom {y, z} follows similarly. 
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Figure shows the numerical solutions of Eq. (fl4j) with the boundary conditions: 
d x ll x (0) = and U x (0) = 1, for several small values of T. The reason for choosing small 
T will be clear later. One can see that the spatial part of quantum probability density 
p x (x) is being trapped by its own self-generated quantum potential U x (x) [23j. Moreover, 
there is a finite distance x = ±x m at which the partial quantum potential is blowing-up 
U x (±x m ) = oo, so that the corresponding partial quantum probability density is vanishing, 



p x (±x m ) = 0. This is a familiar phenomena in nonlinear differential equation 24j, which 
for the case at hand, can be proven as follows. 

Let us define a new variable u x (x) = d x ll x . The nonlinear differential equation of Eq. 
(Il4p then transforms into 

1 AT 

d x u x = —u\ + j^U x . (15) 

The boundary condition translates into u x (0) = d x U x (0) = 0. Further, let us now consider 
the following nonlinear differential equation 

1 2 AT 

d x u x = —u x + —X, (16) 

where X = U x (0); with u x (0) = 0. Since U x (x) > £4(0) = X, then it is obvious that 
\u x {x)\ > \u x (x)\. 




FIG. 1: (a) The x— part of spatial quantum probability density p x (x) (solid line) and its corre- 
sponding partial quantum potential (dashed line) U x (x) for several small values of T obtained by 
solving equation (|14p . We also plot the analytical solution for p Xo (x) at T = 0, assuming that 
U Xo (x) takes the form of a box with infinite wall, (b) u x (x) and u x {x). See text for detail. 
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FIG. 2: (a) The square box spatial support and (b) the half-length of the x— side of the box x m 
plotted against T. 

One can then solve the latter nonlinear differential equation of Eq. (fT6l) analytically to 
have: 

2T 1 

u x {x) = atan(fcr), a = — v2X, b=—v2X. (17) 

It is then clear that at x — ±x m = ±7r/(26), u x is blowing-up, namely u x (±x m ) = ±oo. 
Recalling the fact that |m x (x)| > |w x (x)|, then u x (x) is also blowing-up at points x = ±x m , 
u x (±x m ) = ±oo, where x m < Fig. [Tb. Hence, one can conclude that U x (x) is also 

blowing-up at x = ±x m , U x (±x m ) = oo. It is then safe to say that the x— part of spatial self- 
trapped quantum probability density p x (x) possesses only a finite range of spatial support: 
■M-x = [— x m ,x m ). Finally, the spatial part self-trapped quantum probability density p x (x) 
possesses only finite-size spatial support: .M x = [—x m ,x m ] ® [— y m , y m ] ® [— z m ,z m ], which 
takes the form of a three dimensional square box with sides length 2i m , i — x,y, z. See Fig. 
Hi. 

Let us see what happens if one varies the parameter T. Figure [2b shows the values of 
function of T obtained by numerically solving the differential equation of Eq. ( TH|) 
with fixed boundary conditions: d x ll x (0) = 0, U x (0) = 1. One can see that as we increase 
T, x m decreases, and eventually vanishes for infinite value of T. This shows that p x (x) is 
converging toward a delta function for infinite T. A very interesting fact is seen for the 
opposite limit of vanishing T. One observes that lim-r^o x m (T) = xq, where Xq is finite. 
This fact suggests to us that at T = 0, the spatial quantum probability density and thus its 
corresponding quantum potential are converging toward certain functions: 

hmp x (x;T) = p xo {x), lim U x (x;T) = U xo (x). (18) 
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Let us discuss this latter asymptotic situation in more detail. First, one can see in Fig. 
[T^i that as T decreases, the quantum potential inside the spatial support is getting flatterer 
before becoming infinite at the boundary points, x = ±x m (T). One might then guess that 
at T = 0, the quantum potential is perfectly flat inside the one dimensional box of spatial 
support and is infinite at its boundary points: x = ±x . Guided by this guess, let us 
calculate the profile of the spatial quantum probability density for vanishing value of T. To 
do this, let us denote the assumed positive definite constant value of the quantum potential 
inside the support as U Xc . Recalling the definition of spatial quantum potential U x (x) given 
in Eq. ( 1131) one has 

d% (x) = -^L (x), (19) 

where I XQ = y/p^. The above differential equation must be subjected to the spatial boundary 
condition: p XQ (±x ) = I XQ (±x ) = 0. Solving Eq. (Tl9l one has 

I xo (x) = A xo cos(k xo x), (20) 

where A xo is normalization constant and the wave number k xo is related to the quantum 
potential as: 

k xo = x]W x JW. (21) 

The boundary condition implies 

k Xo x = n/2. (22) 

Figure dh shows that as T decreases toward zero, p x (x;T) obtained by numerically solving 
Eq. (TH|) is indeed converging toward p Xo (x) obtained in Eq. (120]) . This observation thus 
confirms our guess that at T = 0, the x—part spatial quantum potential is flat inside the 
spatial support M. x and is infinite at its boundary points. 

Hence, at T = 0, in total the spatial part of the quantum probability density can be 
written as 

Pxo(x) = IJ Pi (i), (23) 

i=x,y,z 

where pi = A io cos (k io i) , i = x,y,z. Moreover, the support is given by three dimensional 
volume -M Xo = [— Xq,xq] ® [—yo,Uo] ® [~ z 0i z o] at surface of which the spatial quantum 
potential is blowing-up so that the quantum probability density is vanishing. 



S 




FIG. 3: The profile of temporal self-trapped quantum probability density pt(t) (solid line) and its 
corresponding temporal quantum potential Ut(t) (dashed line) for several small values of T. See 
text for detail. 

B. Temporal self-trapping 

Next, let us discuss the temporal part of quantum probability density p t (t) by solving 
the lower differential equation in Eqs. pip . In particular, we are interested to investigate 
the behavior of p t (t) at the limit T — > 0, if it exists. Figure [3] shows the solution with the 
boundary: d t U t (0) = and U t (0) = U tc = — 1. pt(t) and the corresponding U t (t) are plotted 
for several small values of T. One can again see similar phenomena with the spatial part 
that p t {t) is being self-trapped by the corresponding U t (t). One also sees that the support 
of p t (t) is finite given by the interval Ait = [— t a ,t a ] at the boundary points of which the 
temporal quantum potential is blowing-up: U t (±t a ) = oo. 

In Fig. H]we plot the variation of t a against T while keeping U t (0) = U tc fixed. First, in 
contrast to the spatial part, one observes that t a is a monotonically increasing function of T. 
Further, in contrast to the spatial part in which the length of the support is vanishing for 
infinite value of T, the length of the support of temporal quantum probability density p t (t) 
is blowing-up at finite value of T = T p (U tc ). See Fig. HI Given U tc , then for T > T p {U tc ), 
the temporal quantum potential Ut(t) is no more convex everywhere, but is almost periodic 
as depicted in Fig. [51 Hence, p t {t) as defined in Eq. ([6]) is no more normalizable. This case 
is therefore physically irrelevant. 

Yet, again, as in the case of spatial quantum probability density, as one decreases T 
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toward zero, t a is converging toward a finite value to: 

lim t a = t . (24) 
This again shows that at T = 0, U t (t) and p t (t) will converge toward some functions: 

lim U t (t;T) = U t0 (t), hmp t (t;T) = p to (t). (25) 
Below we shall be interested to further study the case of vanishing T. 



t a 



to 



FIG. 4: The variation of half support t a of the self-trapped temporal quantum probability density 
pt(t) against T. See text for detail. 

Proceeding in the same way as for the spatial part, let us calculate pt {t). To do this, 
first one observes in Fig. [3] that as T is approaching zero, U t (t) is getting flatterer inside 
the support before becoming infinite at the boundary points, t = ±t a (T). Again, let us 
guess that at T = 0, the temporal quantum potential is perfectly flat inside the support 
Ait — {—to, to), given by U tc < 0; and is infinite at t = ±t - Recalling the definition of 
U t (t) given in Eq. ffTUl) . one has 

Wto = (26) 
Here I to = J~f%- The above differential equation must be subjected to the boundary condi- 
tion: I to (±to) = 0. Solving Eq. (|26j) . one obtains: 

I to (t) = A to coB(u t). (27) 

Here A to is a normalization constant and the angular frequency uj is related to the quantum 
potential as 

uo = y/(-2<?U tc /h 2 ). (28) 
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FIG. 5: Periodic temporal quantum potential for the case T = 1 with dtUt(0) = and (0) = — 1. 
See text for detail. 

The boundary imposes: 

tu t = 7r/2. (29) 

One finally sees in Fig. [3] that as one decreases T toward zero, p t (t; T) obtained by solving 
the lower differential equation in Eqs. ffTTl) is indeed converging toward Pt {t) given in Eq. 
fl27|) . This again justifies our guess that at T = 0, U t {t) is perfectly flat inside M.t , and is 
infinite at t = ±t . 

III. SPACETIME SOLITON 

Hence, in total, at T = 0, po(q) = p xo (x)p to (t) satisfies the differential equation of Eq. 
(J7J). Notice that the spatio-temporal support of po(q) is composed by M. = -M xo ® M-t - 
Inside A4, the quantum potential is thus flat given by 

U{q) = U Xc + U yc + U Zc + U tc = \ (t?kl - ^) = U , (30) 

where k^_ o = k XQ • k XQ and k XQ = {k XQ , k yo , k ZQ }; and we have employed Eqs. (ED) and (1281) . 
One therefore observes that at T = 0, the quantum force is vanishing inside the spatio- 
temporal support, d a U = 0. 

Now, at r = 0, let us choose the following initial wave function, {po(q), Vo(q)}- Here 
Po(q) = p xo (x)p to (t) and t> a (<z) is a uniform velocity vector field having non- vanishing value 
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only inside the spatio-temporal support A4, that is t> a (<?) = Vq, where Vq is a constant 
four velocity vector. Since at r = the quantum force is vanishing, then initially one has 
dv a /dr = 0. Hence, at infinitesimal lapse of proper time, r = At, the velocity field is kept 
uniform and constant. This in turn will shift the initial quantum probability density in 
spacetime by Aq a = VqAt, while keeping its profile unchanged: p(q; At) = po(q a — VqAt). 
Accordingly, the spatio-temporal support will also be shifted by the same amount: A4at- 
The same thing will happen for the next infinitesimal lapse of proper time and so on and so 
forth. Hence, at finite lapse of proper time r, one concludes that the pair of fields 

{p(q; r), v a (q; r)} = {p (q b - v&r), v a (q b - v&r)}, (31) 

comprises the stationary wave function of the relativistic Madelung fluid dynamics. Here, q 
belongs to the spatio-temporal support at proper time r, denoted by M. T . We have thus a 
spatio-temporally localized wave packet traveling in spacetime: namely a spacetime soliton. 
See Fig. El 
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FIG. 6: Spatio-temporally localized wave packet moving in spacetime. 



A. Mass 

Next, before proceeding to write down the explicit form of the spacetime soliton wave 
function, let us first discuss the physical meaning of quantum potential. First, since the 
quantum potential is constant inside M. r given by Uq, one has U — J dqlfp = Uq. Let us 
proceed to choose a sufficiently large u)q by picking sufficiently large \U tc \ so that U = [7 
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given in Eq. (!30j) is negative, U < 0. This allows us to define a new quantity m as 

C/ = U = ~m 2 (?, (32) 
so that inserting into Eq. (13(1 one obtains 

^> - h 2 kl Q = m 2 c 2 . (33) 
Recalling again the definition of quantum potential of Eq. (JSJ), Eq. (13"2"1) can be rewritten as 

T7i 2 C 2 

□7(g;r) + -^-/(g;r) = 0, 96^, (34) 

where / = J~p. 

A physical interpretation to the above formalism is in order. Eq. (1341) is but the Klein- 
Gordon equation with mass term m. One however should keep in mind that the differential 
equation of Eq. ( 134"|) must be subjected to the boundary condition that I(q; r) is vanishing 
at the surface boundary of the spatio-temporal support, M. T . One can also see that though 
I(q) is not Lorentz invariant, m is. Moreover, since Uq is conserved so is m. 



B. Energy-momentum 

Let us proceed to discuss another conserved quantity: v a = Vc" '■ To do this, let us use 
the conserved invariant quantity m developed in the previous subsection to rescale the affine 
parameter A as follows: 

A = mX. (35) 
One thus has d/dX = md/dX. Using this, Eq. (j2J) becomes 

itM = —Dip. (36) 

Moreover, the pair of equation in Eqs. ([3D translates into 

m^t = -d a U, ^ + d a (pv a )=0, (37) 
dX oX 

where v a = dq a /dX, and is now related to the quantum phase S(q) as 

v a = d a S/m = v a /m, (38) 
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and U is given by 

U = U / m = (39) 
2m I ( J 

Let us then use the proper time r for the affine parameter, r = A, and to avoid complicated 
notation let us rewrite the rescaled four velocity vector v a back as v a . One can then define 
four momentum vector as 

p a = mv a . (40) 

which is conserved. The spatial part of the above quantity gives us the usual definition of 
classical momentum: 

p = mv. (41) 
From the temporal part one can define a scalar quantity E as 

E = cp°= -> mc 2 + -mv 2 , (42) 

y/\ ~ (v/cY 2 

where v is the absolute value of spatial velocity vector, v — ||v||, and the arrow implies that 
v/c is sufficiently small. Hence, E has the dimension of energy. The second term is the 
kinetic energy and the first term is usually dubbed as rest-mass energy, namely the energy 
when the soliton is not moving. 

Next, let us calculate a quantity which is convenient for later discussion defined as follows: 

(F) = J dqr{q){^p)m- (43) 

Since it is conserved, then it is sufficient to calculate its value at r = 0. Putting the wave 
function in polar form ip = I exp(iS /h), one gets 

(P) = [ dq (lo^-ni + ^daSod-So 
J v Z tyi Z i yi 

--I d a I d a S + i^I 2 nS ) . (44) 
m 2m J 

The first term on the right hand side is nothing but the average quantum potential divided 
by mass m: 

[ dqI ^-ni = -[ dqp (q)U (q) = ^. (45) 
Jm 2m m J m m 

Recalling the relation between the four velocity vector field given in Eq. fl38l) . the second 
term on the right hand side is given by 

K = f dqp (p a p a /(2m)) = (46) 
14 



where in the last equality we have used the fact that the velocity vector is uniform inside 
the spatio-temporal support. Further, for the case of soliton where v a = d a So/m is uniform, 
the last term is vanishing due to the vanishing divergence d a v a = — DSo/m = 0. Finally, the 
third term is also vanishing due to the fact that Io(q) = I to (t)I xo (x)I yo (y)I zo (z) is separable 
and each Ii (i) possesses a symmetry J io (z) = I io (—i) so that diI io satisfies 0jij o (— i) = 
—dili (i), i = t,x, y, z. Hence, in total one obtains 

^2, ,2 



-mc 2 . (47) 



( — m 2 c 2 + p a Pa) 2 



2m 

Hence it is given by the negative of the rest-mass energy. 



C. de Broglie's wave function for a single free particle 

Now let us write down the complex- valued spacetime wave function i/j(q; r) = I exp(iS / h) 
corresponding to the spacetime soliton that we have just developed. To do this, one has to 
calculate the quantum phase by integrating d a S = p a to give us 

E 

S(q;r) = -— (ct) +p x x+p y y + p z z + £(T), (48) 

where £ (r) is a function only of r. Hence, at proper time r, one gets 

ip(q; r) = J X0 (x - vr)J to (ci - Er/ (mc)) 

xexp (i(-Et+p x x + v y y + p z z + £(T))/h\ (49) 

Inserting this into the relativistic Schrodinger equation of Eq. (J2]) one can show that £ (r) is 
related to (P) as d^/dr = —(P) which can be integrated to give £(r) = — (P)r = mc 2 r up 
to some constant. Putting this back into Eq. (|4*9|) one finally obtains 

ip(q; t) = Jxo(x - vr)I to (ct - Erj (mc)) 

x exp (i(—Et + p x x + p y y + p z z + mc 2 r)/^j . 

(50) 

One concludes that the phase of the soliton wave function uniquely gives the mass-energy- 
momentum of the particle. 
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One can also see from Eq. (1501) that the soliton wave function possesses internal vibration 
in spacetime whose wave number and angular frequency are given by the four momentum 
vector p a as 

kBi=Pi/fr, i = %,y,z, and uj b = E / h. (51) 

The above relations are nothing but the de Broglie's conjecture in his attempt to explain the 
dual nature of matter as both particle and wave by relating the particle properties (energy- 
momentum) of the matter to the wave properties (angular frequency-wave number) of the 
matter. Yet, in contrast to our approach which is based on linear wave equation, de Broglie 
envisioned such soliton solution to be derived from a nonlinear wave equation. 

Now, let us assume that the internal vibration resonates inside the spatio-temporal sup- 
port. Namely, the wave number of the internal vibration are equal to the integer multiple 
of the wave number of the spatial part of quantum amplitude; and moreover, the angular 
frequency of the internal vibration is equal to the integer multiple of the angular frequency 
of the temporal part of the quantum amplitude: 

k-Bi — Tiiki , i = x,y, z, 
to B = n t u Q . (52) 

where n x ,n y ,n z ,n t are integer. Multiplying both sides of the above equations with Planck 
constant h, one gets 

Pi = hk Bi = nihk io , i = x,y, z, 
E = hu)B = n t fkj . (53) 

The above obtained relations tells us that the energy-momentum are quantized into discrete 
values. 

Further, let us proceed to discuss a special case when n x = n y = n z = n t = 1 to have 

E = hu , p = frk X0 . (54) 

Using the above relation between the energy-momentum of the soliton and the angular 
frequency-wave number of the spacetime quantum amplitude, then Eq. (!33|) translates into: 

E 2 



c 



2 - P 2 = m 2 c 2 . (55) 
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This is the classical energy-momentum relation of special relativity. Since in general one can 
choose k Xo and loq independently of p and E, then the above result opens the possibility of 
the violation of classical energy-momentum relation provided that Eq. (1541) is not satisfied. 

It is then interesting to ask when the resonance of the internal spatio-temporal vibration 
inside the spacetime support occurs. To answer this, it is intuitive to learn from the classical 
wave phenomena. In this case, the resonance phenomena occurs if there is interaction and 
transfer of energy. One therefore might guess that the quantization of energy-momentum of 
Eqs. ( |53|) occurs if there is interaction and transfer of energy-momentum. In other words, 
we expect that the energy-momentum can only be transferred from one matter to the other 
through interaction in discrete quantized value. For the case of light, this statement is but 
the Planck-Einstein conjecture which eventually gave birth to quantum theory. To discuss 
this issue in precise manner, one has to develop a theory describing interaction among 
particles. 

D. Spacetime uncertainties 

One of the important property of the spacetime soliton wave function we developed in 
the previous subsections is that it is broadened with finite width both in space and time 
axes. Namely, any particle should be considered as an extended object both in space and 
time. One is thus suggested to abandon the view of seeing a particle as a dimensionless 
point except for approximation to certain situation. 

It is well-known that the assumption of point particle has led to many formal and physical 
difficulties in physics. It is for example argued as the origin of the divergence of calculations 
in the current version of quantum field theory. Another example is the famous paradox of 
self- interference in double slits experiment [3|]. Assuming a point particle will force one to say 
an ambiguous sentence that the particle is passing through both slits to interfere with itself at 
the screen. Our soliton model of particle developed in this paper might then be considered as 
prospective candidate to address these important issues. It is thus imperative to understand 
how the width of the soliton wave function is related to other directly observable physical 
quantities. 

First, let us discuss the spatial width of the soliton wave function defined by A, = 2i , 
i = x,y, z. Recalling the relation k io i = ir/2, i = x, y, z of Eq. which is coming from 
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the boundary condition at the blowing-up point, one has 

A, = — , i = x,y,z. (56) 

Hence the width in the i— axis is proportional to the inverse of the i— part of the wave number 
of the spatial quantum amplitude. To see its relation with directly observable quantities, it 
is convenient to define the following quantity: 

A A X .A,A 2 

x V( A v) 2 (^) 2 + (a.) 2 (a 2 ) 2 + (A x y(A y y 

IT 7tH 

=• (57) 



'k 2 I h 2 u)i 



x ° y —^2 — m C 

It depends on the angular frequency of the quantum amplitude and also on mass. In the case 
when the internal vibration resonates inside the spatio-temporal support in its first mode, 
the denominator is equal to the absolute value of the momentum of the particle to give: 

x = / = 71 — m- (58) 

ht — m 2 c 2 "P" 

Next, let us discuss the temporal width of the soliton wave function defined similarly by: 
At = 2t . Again, using the relation u t = n/2 of Eq. ff29l) . one gets 

A* = — = — . (59) 

c^h 2 kQ + m 2 c 2 

Similarly, notice that if the internal vibration beats inside the spatio-temporal support in 
its first mode, the denominator is just equal to the energy E to give us 

A, = (60) 



IV. SUPERPOSITION OF MASSES 



Now let us proceed to discuss the implication of the linearity of the Schrodinger equation 
of Eq. Recall that the soliton wave function given in Eq. f)5Up can be interpreted as 
the wave function of a single free particle with a given mass m, energy E and momentum 
P = {P xi Pyi Pz 

}. Since the Schrodinger equation of Eq. (T5]) is linear with respect to the 
wave function, then any superposition of the solutions of the type given in Eq. (!50j) will also 
satisfy the Schrodinger equation. 
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For example, let us assume that ipi{q\ t; {mi, Ei, Pi}), i = 1,2 takes the soliton form given 
in Eq fl50|) with the mass m^, energy- momentum {Ei, pi}, i = 1,2. Then, the following wave 
function which is obtained by superposing the two solitons wave functions: 



ip(q; r) = A x V>i(g; r; {mi, E u pj) 



+A 2 ip 2 (q; t; {m 2 , E 2 , p 2 }) 



(61) 



also satisfies the Schrodinger equation of Eq. ([2]). Here A i: i = 1,2 are two complex- valued 
constants which satisfy: \Ai\ 2 + \A 2 \ 2 = 1. One can thus interpret each soliton term on the 
right hand side as describing a particle with various values of masses, energies and momenta. 

Notice that at period of time when the spatio-temporal support of both solitons wave 
functions are overlapping, one will observe interference in space and time. Hence it is 
impossible to distinguish one soliton from the other. One therefore should consider both 
as a single particle. Now let us assume that pi = — p 2 , namely each soliton is moving 
with direction opposite to the movement of the other. Then at sufficiently large time r, 
the two solitons are no more overlapping in space so that one can already distinguish one 
from the other. Both might however still overlap in time domain. If further one chooses 
Ei/mx ^ E 2 /m 2 , then at sufficiently large proper time, both soliton will also separate in 
time axis. 

Notice also that both solitons move independently from the other. One can thereby 
calculate the total mass, m 12 , to obtain 



This can be shown easily by choosing proper time r so that both solitons are not overlapping 
in spacetime so that one has U = U\+U 2 , where JJ\ is the average quantum potential of soliton 
i = 1,2. Hence, the manifestly covariant Schrodinger equation of Eq. (T5]) can be seen as field 
theory in which a particle of mass m\ 2 can break into two particles of masses m\ and m 2 so 
that mi + m 2 = m\ 2 . Each particle moves independently with momentum pj and possessing 
energy Ei, where i = 1,2. Moreover, since energy and momentum must be conserved then 
one can define the total energy-momentum as Ei 2 = E\ + E 2 and pi 2 = Pi + p 2 - Conversely, 
since the Schrodinger equation of Eq. (T5]) is time reversal, then by reversing the evolution 
of affine parameter, the initially separated two solitons each describing a particle of mass 




(62) 
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m», % = 1, 2, can merge into a single particle with mass = iri\ + vn-i- Needless to say, 
extension into the superposition of more than two solitons are straight forward. 

V. CONCLUSION AND DISCUSSION 

We have thus shown that first, the manifestly covariant mass-less Schrodinger equation of 
Eq. (J2J) possesses a class of non-dispersive soliton solutions with finite-size spatio-temporal 
support. The quantum amplitude inside the support satisfies Klein-Gordon equation with 
finite emergent mass. We finally came to give an interpretation to the soliton solution as 
describing a single particle with finite mass, energy and momentum. Moreover, we showed 
that the soliton solution possesses internal spatio-temporal vibration with angular frequency 
and wave number which are determined by the energy-momentum of the particle as exactly 
conjectured by de Broglie while proposing the solution for the problem of duality of matter 
as both particle and wave. However, in contrast to his envision of nonlinear wave equation 
assuming soliton solution, we showed that our soliton wave function is a solution of a linear 
wave equation. We argued that the fact that the soliton wave function has finite extension 
in spacetime might give the key to the problematic divergence problem encountered in the 
current version of quantum field theory and also to the paradox of particle self-interference 
in double slits experiment. 

We showed that if the internal vibration resonates inside the spatio-temporal support, 
then one can show that the energy and momentum are discretized into packets as assumed by 
Planck and Einstein. We further suggested that this situation appears if there is interaction 
among matters which induces transfer of energy. We also showed that the classical energy- 
momentum relation is recovered when the internal vibration resonates inside the spatio- 
temporal support in its first excited mode. This opens the possibility of violating the energy- 
momentum relation if one moves away from the first resonance mode. 

Next, the linearity of the Schrodinger equation allows for superposition of such solutions 
to further comprise a class of many solitons solutions. Each term of the superposition then 
describes a single particle with a given mass, energy and momentum, moving independently 
of the other. One can thus consider the wave function as a real physical field and regards the 
Schrodinger equation of Eq. (J2J) as a field theory describing multi-particles systems. In this 
theory, things thus live in the ordinary three dimensional space rather than in configuration 
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space. 

There are at least three other interesting issues suggest further exploration. First, we 
have identified soliton as particle by attributing to it mass and energy-momentum. It is then 
natural to ask: what about spin? Second, since the soliton wave function is broadened also in 
time axis, then it might be possible to see interference in the time domain [2^] . Finally third, 
our theory allows us to construct a superposition of two solitons which initially interferes 
each other but then is separated from each other as time goes. It is then interesting to 
employ this idea to re-think the EPR thought experiment. 
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